A dynamic jamming point for shear thickening suspensions 
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We report on rheometry measurements to characterize critical behavior in two model shear thick- 
ening suspensions: cornstarch in water and glass spheres in oil. The slope of the shear thickening 
part of the viscosity curve is found to increase dramatically with packing fraction and diverge at a 
critical packing fraction ^c- The magnitude of the viscosity and the yield stress axe also found to 
have scalings that diverge at (j)c. We observe shear thickening as long as the yield stress is less than 
the stress at the viscosity maximum. Above this point the suspensions transition to purely shear 
thinning. Based on these data we present a dynamic jamming phase diagram for suspensions and 
show that a limiting case of shear thickening corresponds to a jammed state. 



In Newtonian fluids the viscosity does not change with 
apphed shear rate, while non-Newtonian fluids usually 
show a decrease of viscosity when sheared faster, i. e., 
they shear thin. The opposite behavior, shear thicken- 
ing, is less common but can be quite dramatic: beyond 
a certain shear rate the viscosity increases potentially 
by orders of magnitude. This behavior is reversible so 
the stress relaxes when the shear is removed. Reported 
shear thickeningfluids are usually densely packed colloids 
or suspensions [J [3i [ll • Strong shear thickening can be 
observed for example with cornstarch in water. Shear 
thickeningis a concern across a range of industrial pro- 
cesses and is of interest for the ability to absorb 
energy from impacts 0]. 

It has been suggested 0, i, 0, S i, d IH E El, 0, 
[T5I , [TgI . [TtI that shear thickening is related to the phe- 
nomenon of jamming. The concept of a jamming tran- 
sition, however, applies to the limit of vanishing shear 
rate and is exemplifled by the onset of glassy behavior 
with a seemingly diverging viscosity in molecular liquids, 
dense packings of colloids, or macroscopic granular ma- 
terials, leading to the appearance of a yield stress below 
which there is no flow [ll, [l^, . Shear thickening, on 
the other hand, occurs at non-zero shear rate. While a 
yield stress has been measured in some shear thicken- 
ing fluids [1], no linkage has been established between 
such yield stress and the observed shear thickening be- 
havior. Dense shear thickening fluids have been reported 
to exhibit seemingly discontinuous jurnps in stress with 
increasing shear rate 0, Q [13, HH [13, [ll [13, HI] . How- 
ever, this discontinuity has not yet been characterized 
quantitatively. Phenomenological shear thickening mod- 
els have shown stress distributions that are similar to 
those of force chains in jammed systems [1, [l3|. How- 
ever, the mild shear thickening behavior found in simula- 
tions [13, 126,, ^8] so far cannot reproduce the dramatic 
increases in viscosity with shear rate observed in many 
experiments. As a result, the connection between shear 
thickening and the onset of jamming has been qualitative 
and many details of the relationship remain unresolved. 

Here we characterize the stress-shear rate discontinu- 
ity and the relationship to the yield stress through rheo- 
logical measurements. This includes the first systematic 
characterization of shear thickening rheology as a func- 



tion of packing fraction for suspensions of non-Brownian 
particles. We find that discontinuous shear thickening is 
a limiting behavior which is approached at a critical pack- 
ing fraction where the onset shear rate of shear thickening 
approaches zero and the yield stress jumps dramatically. 
In other words, the limiting case of discontinuous shear 
thickening corresponds to a jammed state. We then de- 
velop a phase diagram to delineate the shear thickening 
and jammed regimes. 

We present below experimental results from two rather 
different hard particle suspensions, cornstarch in water 
and glass spheres in mineral oil, to demonstrate the gen- 
erality of our results. Cornstarch (average particle di- 
ameter of 14 /xm) was chosen as a prototypical shear 
thickener, while glass spheres (88-125 /im diameter) have 
the advantages of better defined particle properties and 
better studied packing properties. Measurements were 
performed in a rheometer using Couette (inner rotat- 
ing cylinder of diameter 26.6 mm, gap of 1.13 mm) or 
parallel plate (25 mm diameter rotating top plate) ge- 
ometries. All data presented were taken with a constant 
applied torque converted to a shear stress r (based on a 
Newtonian fiow profile). The shear rate 7 is defined by 
the measured rotation velocity over gap size. We define 
viscosity by ry = r/j. Samples were pre-sheared before 
experiments and measurements were performed with in- 
creasing as well as decreasing stress ramps to check for 
repeatability. Reported packing fractions are based on 
measured particle and fluid masses mixed together before 
shearing. 

The starch particles were suspended in water, density- 
matched to 1.59 g/mL by dissolving CsCl. From mea- 
surements with optical tweezers we found no significant 
particle interactions in the stress range of the experi- 
ments. We measured packing fractions using the weight 
of starch at ambient conditions of 23° C and 42% humid- 
ity which includes some water weight. Glass spheres with 
hydrophobic coating were used to optimize dispersion in 
mineral oil (viscosity 58 mPa-s) to minimize the yield 
stress. Since glass spheres are denser than mineral oil, 
measurements were performed in the parallel plate geom- 
etry with a gap size of 0.5 mm. This minimized pressure 
on the packing due to its own weight. Measurements re- 
ported here are done with no oil extending outside the 
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FIG. 1: Shear stress r vs. shear rate 7 for packing fractions 
(j), as indicated (color online), (a) Glass spheres in mineral 
oil. (b) Cornstarch in water. Dashed lines: slope 1 indicating 
fixed viscosity. 



parallel plates so that the particles are confined to the 
space between the plates by surface tension. 

Figure [1] gives the evolution of stress vs. shear rate 
curves as the packing fraction is increased. On a log-log 
plot as shown, a slope of 1 corresponds to Newtonian 
flow, indicated by dashed lines for reference, a slope be- 
tween zero and unity corresponds to shear thinning, while 
a slope greater than unity signals shear thickening. The 
overall steepness of the traces within the shear thicken- 
ing region is seen to increase with (f> and to approach a 
vertical line where shear thickening becomes discontinu- 
ous. Another feature is that shear thickening is associ- 
ated with a certain intermediate stress range that varies 
little with packing fraction. For stresses either larger 
than the upper limit of the shear thickening region or 
smaller than the shear thickening onset, shear thinning 
behavior is observed. At lower <j) the slopes gradually 
approach 1 at all stress ranges so there is a gradual tran- 
sition to Newtonian flow. The behavior described above 
is similar to what has been found by careful measure- 
ments in shear thickening colloids [251 . The similarity 
is notable because it is usually assumed that Brownian 
motion and electrostatics are important factors in shear 
thickening ^ which are insignificant for our larger parti- 
cles. At sufficiently large 0, the curves exhibit a non-zero 
stress value in the limit of zero shear rate, i. e., a yield 
stress. Given our stress resolution around 10~^ Pa in 
our rheometer, this is most clearly seen in the cornstarch 
data where the yield stress is larger. The yield stress is 
seen to encroach on the shear thickening stress range at 
high packing fractions above which there is only shear 
thinning and no shear thickening. 

To quantify the stress-strain relationship in the shear 
thickening regions, we fit the traces locally to a power 
law rj oc T^~'^ which is equivalent to r oc 7^/'^ but can 
be fit because its slope does not diverge. The param- 
eter e depends on packing fraction and corresponds to 
the inverse slope of the traces in Fig. [TJ Newtonian flow 
corresponds to e = 1 and a stress divergence corresponds 
to e = 0. In Figs. 2a&b, the e values plotted are from 
fits around the steepest portions of the stress-shear rate 




FIG. 2: Evolution of viscosity and yield stress with packing 
fraction <f) for glass spheres in mineral oil (a,c,e) and corn- 
starch in water (b,d,f). (a,b) Inverse logarithmic slope e of 
stress-shear rate traces in the shear thickening regime, de- 
fined by T (X 7^^'^. Insets show micrographs of the particles. 
(c,d) Viscosity scale fj defined as geometric mean of the vis- 
cosity in the shear thickening region. Insets in (c,d): log-log 
plots of the same data relative to ^c. (e,f) Yield stress Ty. 
Solid triangles: increasing stress measurements . Open tri- 
angles: decreasing stress measurements. Dotted lines: the 
plateau yield stress Ty^max- Inset in (e): detail view of the 
region close to (f>c. All panels: Solid lines are power law fits 
explained in text, dashed lines are resulting values of (j)c. 



traces. For both starch particles and glass spheres, e ap- 
proaches zero at a critical packing fraction (pc where the 
slope of the viscosity curve becomes divergent. Previous 
qualitative descriptions of shear thickening have made a 
distinction between continuous and discontinuous shear 
thickening regimes, implying the possibility of e = over 
a range of cf). However, the fact that e only approaches 
zero at (pc suggests that discontinuous shear thickening 
is better thought of as a limiting behavior of shear thick- 
ening. The value of 0c is obtained from a power law fit 
of e oc (0c — 0)" to the data in Figs, d^&b which gives 
0c = 0.478 ± 0.003 and n = 0.5 ± 0.2 for starch and 
0c = 0.564 ± 0.004 and n = 0.5 ± 0.2 for glass with rel- 
ative statistical uncertainties. Power laws shown are fit 
from 0c down to the smallest which is consistent with 
a power law fit. When the lower end of the fit range 
of was increased, fit values of 0c remained consistent 
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FIG. 3: Dynamic jamming phase diagram for shear thick- 
ening fluids. Data shown is for glass spheres in mineral oiL 
The sohd line corresponds to the yield stress in Fig. [2^, while 
the dotted lines indicate transitions between shear thinning 
and thickening. Solid circle: dynamic jamming point at (pc- 
Open circles: yield stress. Up-pointing triangles: onset of 
shear thickening. Down-pointing triangles: shear thickening 
maximum. 



within quoted uncertainties, however fit values for expo- 
nents varied with the fit range, so we do not claim to 
have measured limiting scaling exponents with certainty. 

As a second indicator of a transition near (pc, we in- 
vestigate the evolution of magnitude of the viscosity. To 
this end, we define a characteristic magnitude fj by the 
geometric mean of the viscosity over a fixed stress range 
in the shear thickening region. Figures 2c&d show that 
fj appears to diverge very close to the point where e goes 
to zero. A diverging power law fj oc (0c — f/*) " is fit to 
the data in Figs. [It&d which gives 0^ = 0.488 ± 0.004 
and n = 3.1 ± 0.3 for starch and (j>c = 0.576 ± 0.004 and 
n = 1.2 ± 0.2 for glass. While the viscosity varies with 
stress the fit value of 0c is independent of the fixed stress 
range chosen for the fit in the shear thickening region. 
This divergent scaling is independent of the yield stress 
because it is in a much higher stress range, and is inde- 
pendent of the diverging slope because that has a mild 
effect on a fixed stress range. 

The yield stress Ty for different packing fractions is 
shown in Figs.[2^&f and is found to increase precipitously 
as the same value 0c is approached. Around 0c, the yield 
stress plateaus at a value Ty^max and does not change 
significantly at higher packing fractions. A power law 
Ty OC (0c — 0)~" is fit to the data in Figs. dJj&f below the 
plateau to obtain 0c = 0.487 ± 0.003 and n = 2.5 ± 0.3 
for starch and 0c — 0.578 ± 0.004 for glass spheres (n 
could not be fit with any certainty because the jump in 
yield stress is so dramatic) . Comparison of the packing 
fraction dependence of e, fj and Ty in Fig. 2 shows that, 
statistically consistent within experimental uncertainties, 
a single value 0c captures the behavior of all three pa- 
rameters for each of the two suspensions 

From the data in Figs. [1^42] we can assemble a phase 
diagram that delineates shear thinning, shear thicken- 
ing and jammed regions in a parameter space given by 
the applied stress Ty and the packing fraction 0. Fig- 



ure [3] shows this for the glass spheres. The jammed re- 
gion here is defined as the portion in the diagram be- 
low the yield stress. As in the usual jamming phase 
diagram fl^, exceeding the yield stress leads to shear 
thinning flow. For < 0c, a new feature in Fig. [2] are 
two boundaries that extend out from the jammed region 
to smaller packing fractions. These boundaries separate 
shear thickening from shear thinning regimes. Increasing 
stress at fixed will take a sample from shear thinning to 
thickening and back to thinning again, as seen from the 
slope changes of T{j) in Fig.[TJ As 0c is approached from 
below, however, the increasing yield stress pushes both 
boundaries upward, forcing them to approach the yield 
stress line that delineates the jammed region. Once the 
yield stress dominates the total stress in the system the 
suspension no longer exhibits shear thickening behavior 
but only jamming and shear thinning. Because of the di- 
vergent scaling of the yield stress, this transition occurs 
at a packing fraction near to but less than 0c. Since the 
onset stress for shear thickening (up-pointing triangles in 
Fig. 3) lies above the yield line, the shear thickening and 
jamming regions are separated by a thin wedge of a shear 
thinning regime. The convergence of the yield stress and 
the boundaries between shear thickening and shear thin- 
ning suggests that 0c is a singular point. We name this 
the dynamic jamming point because it is analogous to a 
static jammed state with a yield stress in the sense that 
the slope divergence implies the stress increases without 
an increase in shear rate, i. e., exhibits a jump in r^j). 
Furthermore, since the viscosity magnitude diverges in 
(j)c — <t> and the onset shear stress of the shear thickening 
regime is relatively fixed except for the influence of the 
yield stress the shear rate required for shear thickening 
extrapolates to zero at 0c. Thus, the discontinuous shear 
thickening limit corresponds to a jammed state. We note, 
however, that the approach to zero onset shear rate can- 
not be measured all the way up to 0c because the diverg- 
ing yield stress suppresses shear thickening before 0c is 
actually reached. This extrapolation is made without ref- 
erence to the yield stress measurements; the observation 
that this dynamic jamming point occurs on the jamming 
phase boundary confirms the correspondence. 

The phase diagram in Fig. [3] is valid for confined sus- 
pensions of non-attracting particles. While details dif- 
fer for cornstarch and glass sphere systems, both show 
qualitatively similar behavior so the overall delineation of 
jammed, shear thickening and shear thinning regions is 
robust. The data shown were taken with smooth rheome- 
ter plates, but we also did experiments with roughened 
plates and saw no qualitative differences. However, there 
is an important factor that can change the phase dia- 
gram as drawn in Fig. [31 First, with attractive interac- 
tions particles can chain across the system to p roduce a 
large yield stress at low packing fractions [23] ■ Such a 
yield stress can overwhelm any stress increase associated 
with shear thickening so that shear thickening is not ob- 
served, similar to the situation above 0c in Fig. [3| The 
boundary conditions are important as well. In the par- 
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allel plate setup, when the experiment was done instead 
with extra fluid extending outside the plates the shear 
thickening curves were less steep, and a divergent scal- 
ing of the slope of the viscosity curve was not achieved. 
This suggests that confinement (in this case by surface 
tension) is necessary to observe the approach to discon- 
tinuous shear thickening. 

Finally, we turn to the significance of the measured 
value for <pc- Given the large uncertainty in determining 
an absolute packing fraction for starch samples due to 
water absorption from the atmosphere, due to the fact 
that dry powders generally pack less tightly than larger 
particles because of van der Waals forces, and due to the 
non-spherical polydisperse shapes of the starch particles 
(Fig. 2b), we focus here on the glass spheres. We find 
near identical values for (pc when the fluid is replaced by 
silicone oil or water, or when we double the particle size. 
Thus, 4>c is not determined by fluid properties or particle 
diameter but rather seems to be related to an issue of ge- 
ometric packing. Our average value of 4>c = 0.573 ±0.007 
(including an absolute error of 0.005) for glass spheres is 
close to the value of random loose packing at 0.56 for 
glass spheres, corresponding to the packing fraction par- 
ticles will reach when settled slowly in the limit of a den- 
sity matched fluid [33|. This coincides with the onset of 
dilation, which is the lowest packing fraction where the 



packing must expand when sheared[30|. The value 0.56 
is for large samples with density matching, and a positive 
correction of around 1% would be expected both because 
of the finite size of our sample and the weight of the pack- 
ing [sij . It is notable that shear thickening has long been 
associated with dilation [il 113, HH 113] , mostly through 
normal force measurements, but a critical packing frac- 
tion has not been identified by such measurements. 

The existence of divergent scalings at a critical packing 
fraction is provocatively reminiscent of a 2nd order phase 
transition. While the jamming point has been shown 
to have similarities to a critical point with some non- 
universal critical exponents [2^1 , and some scalings com- 
patible with a 2nd order phase transition [32] , we are not 
aware of any such model that includes shear thickening. 
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